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ABSTRACT 

We investigate the static solutions of Callan and Maldecena and Gibbons to lowest order 
Dirac-Born-Infeld theory. Among them are charged wormhole solutions connecting branes 
to anti-branes. It is seen that there are no such solutions when the separation between the 
brane and anti-brane is smaller than some minimum value. The minimum distance coincides 
with the energy minimum, and depends monotonically on the charge. Making the charge 
sufficiently large, such that the minimum separation is much bigger than ^fa' , may suppress 
known quantum processes leading to decay of the brane-anti-brane system. For this to be 
possible the zeroth order wormhole solutions should be reasonable approximations of solutions 
in the full D— brane theory. With this in mind we address the question of whether the zeroth 
order solutions are stable under inclusion of higher order corrections to the Dirac-Born-Infeld 
action. 
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1 Introduction 



The Born-Infeld nonlinear description of electrodynamics[l] and its subsequent generalization 
to membranes a la Dirac[2] is of current interest due to its role as an effective theory for 
Dp— branes. The associated Dirac-Born-Infeld (DBI) action appears at lowest order in the 
derivative expansion for the effective Dp-brane action. [3] , [4] The original Born-Infeld theory has 

a charged static solution, which was generalized by Callan and Maldecena[5] and Gibbons [6] 
to families of static solutions on the brane. The families are associated with orbits of the 
5*0(1, 1) group. The Lagrangian is invariant under 50(1, 1) and can be used to label the 
orbits. One such orbit contains the solution of Born and Infeld. Another is a BPS solution 
representing a fundamental string attached to the brane. Finally, there is a third family of 
solutions corresponding to wormholes which connect the brane to an anti-brane a finite distance 
away. Here we show that there are no charged wormhole solutions having a separation between 
the p— brane and anti—p— brane smaller than some minimum value. The minimum separation 



distance goes like |Q|~, Q being the C/(l) charge. At minimum separation, the self-energy 
is also a minimum, where there appears a cusp singularity in the plot of the energy versus 
separation distance. In the quantum analysis of the brane-anti-brane system an instability is 
known to occur at distance scales of order ^/a' due to excitation of tachyonic modes. [7] Then 
for sufficiently large charge, i.e. 



such quantum processes may be suppressed, and it is possible that charged wormholes can 
stabilize the brane-anti-brane system. 

For the above scenario to be correct, however, classical stability of the wormhole solutions 
should be checked. This means a) enlarging to time dependent solutions, and investigating 
whether solutions are stable with respect to perturbations about the static solution. It also 
means b) checking whether the solutions to the zeroth order effective theory are a reasonable 
approximation of solutions to the full effective Dp— brane action. Here we shall only consider 
b). One signal that solutions may be unstable in the sence b) is the presence of singularities 
in the field strength, where the derivative expansion cannot be trusted. Such a singularity 
is present for the original Born-Infeld solution, in that case at a single point, and for the 
entire orbit of solutions connected to the Born-Infeld solution. Despite the singularity, these 
solutions are associated with a finite self-energy. For the BPS case the singularity occurs 
an infinite distance away from the brane, and appears harmless. The wormhole-type solutions 
were originally constructed by joining together two local solutions, obtained in the static gauge, 
at the minimum circumference of the wormhole. [5] , [6] A singularity in the field strength occurs 
along the throat - precisely at the minimum circumference. The singularity in the field strength 
is a coordinate singularity, which can be removed by going to another gauge. Nevertheless, it 
is a signal that higher order corrections may not be negligible. 

To check stability in the sence b), we will rely on recent computations of the derivative 
corrections to Born-Infeld theory. The first order corrections to the action were obtained 




(1.1) 
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separately by Wyllard [8] and Das, Mukhi and Suryanarayana[9]. [Higher order corrections 
seem currently out of reach.] Using their results we carried out a stability check previously 
for the case of the original Born-Infeld solution. [10] There we argued that the original Born- 
Infeld solution is unstable under inclusion of these first order corrections. More specifically, we 
numerically obtained corrections to the zeroth order Born-Infeld solution, but found that they 
give an infinitely large correction to the Lagrangian. We give a simpler proof of the result here. 
Because the Lagrangian is 50(1, 1) invariant the result applies to the entire orbit of solutions 
connected to the Born-Infeld solution. Concerning the BPS solution, it is known that the 
such a solution is stable to all orders in the derivative expansion. [11] We verify that this is 
consistent with the first order results of [8] and [9] . We find that the stability analysis for the 
wormhole solutions leads to the same results obtained for the Born-Infeld solution. Namely, 
corrections to the zeroth order solution lead to an infinitely large correction of the Lagrangian. 
In this case, we need to rely on numerical computation for the result. 

In section 2 we give analytic expressions for the three families of solutions along with their 
self-energies. Here we show that the charged wormhole solutions have a minimum length. The 
question of stability of the zeroth order solutions is addressed in section 3. In appendix A we 
write down the wormhole solution in another gauge, where the field strength is singularity-free. 
In fact, it is a constant in that gauge. In appendix B we use the results of [8] and [9] to obtain 
the first order corrections to the field equations for the BPS case, and show that the answer 
agrees with [11]. 

2 Zeroth Order Solutions 
2.1 Dirac-Born-Infeld theory 

We consider the p-dimensional brane embedded in a ten dimensional space-time with flat metric 
[vab] =diag(— 1, 1, 1), A, B, .. = 0, 1, ...,9. We denote the brane coordinates by X^. They 
are functions of p -|- 1 parameters jjl, u, ... = 0, 1, ...,p. Additional degrees of freedom on the 
brane are ^7(1) potentials Afj,{$,). The DBI action is written in terms of the {p + 1) x {p + 1) 
matrix 

V = ^AB di^X^iO d.X^'iO + 2W F^,{0 , (2.1) 

where df^ = The first term is the induced metric on the brane, while F^i, = d^Av — duAjj, 
is the U{1) field strength. We assume the two-form contribution is absent Bab = 0. The DBI 
action is[l],[2] 

^Sli = ^ [ d^^'C 4°^ , = 1 - J-det[h^u] , (2.2) 

ys J 

where Tp is the tension, which is expressable in terms of a' according to 

Tp = (47r2a')-(^+^)/^ , (2.3) 
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and Qg is the string couphng. ^j^^j is invariant under diffconiorphisms on the branc 
^'^{0 5 ^(1) g^'Ugs transformations A^{^) + 9i^A{^), as weh as ten-dimensional 

Poincare transformations. Variations in X'^{^) and -4^(0 lead to the equations of motion 

5^|V^d^(/i'"' + /i"'^)r7AB5,X^| = 

di,!^V-deth{h'"' -h"")^ = 0, (2.4) 

respectively, where h^^h^p = Sp. 

The known families of spherically symmetric static solutions [5], [6] can be classified in terms 
of 50(1, 1) orbits (we do this below), and they describe different topologies embedded in the flat 
ten-dimensional background. For one family of solutions, containing the original Born-Infeld 
solution, a time slice is W minus a point. Wc could therefore describe it with the introduction 
of a delta function source to the rigiit-hand-side of (2.4). Another family corresponds to a 
brane and anti-brane connected by a worniholc. In that case one can patch together local 
solutions to (2.4). The families of solutions arc written in terms of two integration constants Q 
and C, Q being the electric charge. Locally, all solutions can be expressed in the static gauge, 
where one identifies with the first 1 brane coordinates X^^, n = 0,1, ...,p. The remaining 
X", a = p + l,p + 2, 9, then denote normal coordinates, and (2.1) becomes 

V = V^.u + + 27ra' F^, , (2.5) 

The static solutions of [5], [6] are for a radial electric field with a single transverse mode excited. 
Choose the nonvanishing degrees of freedom to be ^o(^) and Xp^i{r), where r is the radial 
coordinate on the brane. Since the metric is diagonal, the resulting matrix h is diagonal 
except for the 2x2 submatrix with corresponding indices /j, and v equal to and r. That 2x2 
submatrix and its inverse are given by 

-1 -fir) \ 1 f-l-girf -/(r)\ 

fir) l + 5(r)V l + 9{r)'-f{r)'\ f{r) 1 J' 

respectively, where f{r) = 2-iTa' drAo{r) and g{r) = drXp^i{r). Substituting into the equations 
of motion (2.4) gives 

«r ^-^^ . '"'"'^'^ (27) 

The solutions for /(r) and g{r) are 

/(r) gir) 1 



Q C ^Q2_c^^r'^p-2 
The integration constants Q and C have units of [length]^" ^. 



(2.8) 



For the configurations (2.8) the Lagrangian and equations of motion are invariant under 
the SO{l,\) transformation 

/(r) \ / cosh V sinh v\ff{r)\ 
g(r) J \ sinh u cosh v J \ g(r) J 
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The integration constants transform in the same way 

Q\ / coshi/ sinhzy\ / Q\ 

C J \smhu cosh u J \C J ^ ' ' 

There are then three kinds of orbits: i) \Q\ > \C\, ii) \Q\ = \C\ and in) \Q\ < \C\. i) is 
connected to the original Born-Infeld solution, ii) is the BPS solution and Hi) is associated 
with wormhole solutions. The orbits can be classified by their corresponding value for the 
spatial integral of the Lagrangian density C^^^j 



Lg^^.(r) = / dr , / drrP-^ \ 



(2.11) 



where fip-i is the volume of a unit p — 1 sphere. We removed a hole around the origin of 
radius r in the integration for the first integral, in order to accommodate the different cases, 
as their domains differ. The second integral is the vacuum subtraction. The first integral can 
be expressed in terms of the one for Xp^i(r), yielding 

L%ir) = ^!^r^r^P-2 + Q2_c^ + (^2 _ q2)^£^| ^2.12) 

For cases i), ii) and Hi) we find that the result (after setting r to its appropriate value) is 
positive, zero and negative, respectively. 



2.2 Analytic solutions 

The right hand side of (2.8) can be integrated to obtain analytic expressions for the potential 

2p — 2 

.Ao(r) and the transverse displacement Xp^i{r). For this expand in powers of Qi_Qi and 
integrate term by term. The result for the indefinite integral can be expressed in terms of a 
hypergeometric function: 



- C2 V2'2p-2'2p-2' Q2_c'2 



2p — 2 

for \ q-2_(-;2 I < 1- Now set the limits of integration to be r and oo, with the assumption that 
the potentials vanish at the latter. To evaluate (2.13) for these limits we analytic continue 
F(a, b; c; z) using[12] 

Fia,b;c;z) = l^^^^^f ~ "| (-z)-° F{a, 1 - c + a;l - b + a;-) 

r(o)r(c — a) z 

where |arg(— 2;)| < tt. So (2.13) can be rewritten as 

1 f(- y-2 . 3p-4 Q2_c,2 x ^ (Q2_c,2)^ ^/ 2p_i \ /p_2 



(2-p)rP-2 \^2'2p-2'2p-2' r2p-2 J ^ \2p-2J \2p-2 
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since F{a,0;c;z) = 1. For p > 2, only the last term survives when evaluating at oo, which 
then get subtracted out after evaluating between r and cxd. For the transverse displacement 
one gets 



Xp+i{r) 



-C 



{p - 2)rP- 



^ 1 p-2 3p-A Q^-C^ 



2 2p- 2' 2p- 2' 



for p > 2 



(2.15) 



Next we examine the result for the three different orbits: 



Case i) \Q\ > \C\. The limit C — > where the transverse mode is not excited gives the 
original Born-Infeld solution[l], while \Q\ > |C| > yields a deformation of the Born-Infeld 
solution where a spike protrudes from the brane. We plot below the function Xp_|-i(r) for p = 3 
on a two dimensional spatial slice: 




fig.l p = 3 \Q\ = 1 \C\ = .8 
^Prom (2.15) [and (2.14)] the maximum size of the spike is the absolute value of 



(2 -p)Vvr 

For the integral of the Lagrangian density [C.f. (2.12)] one gets 



for p > 2 



^.-1 r(|^)r(^) (Q2_^2)^ ^ 0^ 

p[p - 2)V7r 



(2.16) 



(2.17) 



Case a) \Q\ = \C\. One arrives at the BPS solution for this case, where (2.8) reduce to the 
Coulumb solutions /(r) = g{r) = Since F{a,b;c;0) = 1, (2.15) reduces to 



-Q 



{p - 2) ' 



for p > 2 



(2.18) 



and the spike becomes infinitely long, 
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fig.2 p = 3 \Q\ = \C\ = l 

representing a fundamental string attached to the brane.[5] In this case the integral of the 
Lagrangian density L^^j-(O) goes to zero. These solutions are BPS because they preserve half 
of the supersymmetries of the ground state solution. Supersymmetries are present when the 
matrix 

[r^rf+i] + 2™' [r^n (2.19) 

is degenerate. are F matrices for the ten-dimensional background space, [F'^jF'^] = 2ri^^ . 
To see that this holds when \Q\ = \C\, one observes that (2.19) is proportional to (QF*^ + 
CFP+^)F'", whose square is (C^ - Q^)]!. (2.19) is then nilpotent when \Q\ = \C\. 

Case Hi) \Q\ < \C\. Here one gets a finite diameter tube with a minimum radius vq = 




fig.3 p = 3 \Q\ = l |C| = 1.2 



7 



Both g and /, and consequently also the electric field, are singular at r = rg. Nevertheless, Aq 
and Xp^i are not. From the latter the tube has a finite length. After expressing C in terms 
of ro and Q, it is 



X,+i{ro) = ^ ^ , for p > 2 , (2.20) 
(2-p)r5 il2F^J 

The domain of integration for the integral of the Lagrangian density [C.f. (2.12)] now goes 
from ro to DO. One gets 

^DBiiro) = --J- — ^ ^, % . < , forp > 2 (2.21) 



p(p-2) r( 



The static gauge breaks down at r = ro, and so the above solution is only local. A global 
solution was proposed by gluing this one to the analogous solution on an anti-brane. [5] , [6] 
The global solution then represents a wormhole connecting the brane with an anti-brane a 
distance of 2|Xp-|_i(ro)j away, with a throat of minimum radius tq . The gluing of the two 
local solutions to form a wormhole occurs at r = tq, precisely where there is a singularity in 
the electric field, which might be a matter of concern. On the other hand, the electric field 
singularity is a coordinate singularity, which is easily seen by transforming to another gauge. 
Take for example the gauge where the r— coordinate of the static gauge is replaced by 2; = Xp+i. 
In the new gauge the solution is described by the inverse, call it R[z)^ of the function Xp+i(r). 
Now the electrostatic field is in the z— direction, and there are coordinate singularities at the 
location of the brane (and anti-brane). We denote the electric field (x 27ra') in the new gauge 
by E{z). It can be computed locally by performing a coordinate transformation from the static 
gauge, 

Ei.) = _ = _ m^)) (2.22) 

Substituting in the solution for / and g given in (2.8) gives a constant electric field 

E{z) = Q/C (2.23) 

So in this coordinate frame there are no singularities in the electric field (for C 7^ 0). In ap- 
pendix A we write down the action in this gauge and show that (2.23) solves the corresponding 
equations of motion. 

^Prom (2.20) it follows that there is a minimum separation distance between the brane and 
anti-brane for a fixed Q (and p> 2). It is equal to 

nnn 2|X,„(.„)| = ^ |Q1 A , (2.24) 

(p-2)2p-2 ^\2p-2> 

and occurs when C and Q are constrained by 

= (p - 1)Q2 (2.25) 
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Below we plot the separation distance versus throat size for a fixed Q when p = 3: 




1 2 3 4 5 6 7 

fig.4 p = 3 Q = l 2|Xp+i(ro)| vs. ro 

2.3 Self-Energy 

Concerning the energy, one can apply the canonical formalism starting from the Lagrangian 
in (2.2). For this again assume the static gauge and hence (2.5). The Hamiltonian density is 

= P"Xa + - C%^Bi ' (2-26) 

where the dot denotes a time derivative and 

W = -TraV- det h {h''^ - /i°^) , (2.27) 

are the momenta conjugate to and A^, respectively. As usual, the momentum conjugate 
to Ao is constrained to be zero. After integrating by parts 

'^DBi = P"^o. + + V-deth - 1 - diUUo , (2.28) 

where i = 1,2, ...,p. The coefficient of .Ao gives the Gauss law constraint. The remaining terms 
are equal to 

-V-dcth - 1 (2.29) 

Although the Lagrangian density is invariant under the 50(1,1) symmetry (2.9), the Hamil- 
tonian density is not. Then unlike the integral of the Lagrangian density, the integral of the 
energy density will not be constant along the orbits of S0{1, 1). 

The integral of (2.29) gives the self-energy of the DBI solutions (2.8). After removing a 
hole around the origin of radius r in the integration domain one gets 

m = Sir) , 8{r) = r dr^^^^^l= - H drr^-' , (2.30) 

where the factor Tp/gg comes from (2.2). In the second integral we subtract off the total 
vacuum energy of the brane. Note that we must restrict the lower limit in the first integral to 
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be greater than or equal to tq for the case \C\ > \Q\. The result can be expressed in terms of 

JC... I 1 r , 

(2.31) 



p C p 

This gives a positive answer for the self-energy since Xp+i(r)/C is negative for the solutions, 
while the second term vanishes after evaluating at the minimum value of r (0 for |Q| > |C| 
and ro for \C\ > \Q\). 

For case i) \Q\ > \C\, one gets the total self-energy of the solution by setting r in (2.31) to 
zero, which yields 



U{0) 



T,%-, r(|^)r(2^) (p_i)Q2 + ^2 



2_ ; 

2p-2 



9s p{p-2)^ (Q2_(72)i 

For a fixed Q it goes monotonically from the Born-Infeld value 



for p > 2 



(2.32) 



orn— Infeld 



p-i ^/3p-4 



^ \Q\— 



for p > 2 (2.33) 



gs pip -2)^ \2p-2j \2p-2^ 

corresponding to C = 0, to infinity in the BPS limit, corresponding to |C| — \Q\. We plot 
below f for Q = 1 and p = 3: 



0.2 0.4 0.6 0.8 1 

f ig.5 p = 3 Q = l S vs. C {C <Q) 

For case ii) \Q\ = \C\ the total self-energy is infinite. At large distances the energy 

per unit length of the infinite string solution is constant. From (2.31) it follows that 



\Q\, as \Xp+i\^oo (2.34) 



For case in) \Q\ < \C\, one gets the total sclf-cncrgy by setting r in (2.31) equal to ro: 

For a fixed Q, the minimum energy configuration occurs for v'q' ^ = (p — 2)Q^, corresponding 
to the minimum separation distance (2.24) between the brane and anti-brane. The minimum 
value for U{ro) is 

= 5^ ^^^-^ £|l|i' IQ\^ , tor p > 2 (2.36) 

9s p(p_2)2p-2 ^{2p-2) 
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If (5 = the minimum energy configuration occurs wlien tlie brane and anti-brane coincide. 
For (5 7^ and a separation distance greater than the minimum value (2.24), there are two 
possible solutions with different throat sizes. The one with a smaller throat is energetically 
favored. Call Wo(^p+i) and Ui{Xp^i) the energy of the thin and fat wormholes, respectively. 
For a large separation distance, 

i^oi^p+i) — ^ P P ^ \QXp^i\ 
9s 

Upon plotting the energy versus the separation distance one gets a double-valued function, 
with a cusp at the minimum separation, as is illustrated below for Q = 1 and p = 3: 




4 5 6 7 8 



fig.6 p = 3 Q = l S vs. 2\Xp+i{ro)\ 

The minimum energy solution for fixed Q in case i) \Q\ > \C\ was the original Born-Infeld 
solution, while in case in) \Q\ < \C\ it corresponded to (2.36). In both cases the energy goes 
like IQI^^i. Assuming charge conservation, such solutions are energetically unstable under 
fission into far separated solutions with total charge equal to Q. It was however pointed out in 
[6] that fission may not be realized at the classical level for singular field configurations, and 
the above solutions are of this type. Assuming fission does occur, either classically or quantum 
mechanically the minimum energy configuration should be an ensemble of far separated worm- 
holes in case in) or Born-Infeld solutions in case i) with the fundamental charge. In comparing 
i) with Hi), the ratio Wmin/WBorn-infeid is less than one for p > 4, while it is greater than one 
for p = 3. Thus for p > 4, it is energetically more favorable for wormholes to develop between 
a charged brane and equally charged anti-brane than for Born-Infeld configurations to develop 
on the brane and anti-brane. The opposite is true for p = 3. 

2.4 Thermodynamic considerations 

Here we make a side remark concerning the thermodynamics of wormholes. Once again, for 
case Hi) when the energy is greater than the minimum, two types of wormholes with different 
thickness may be present. Say they are in a heat bath with temperature T and call po and pi 
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the density of the thin and fat wormholes, respectively. If one assumes they are in dissipative 
and thermal equilibrium, then the ratio of their densities at a temperature T is given by 



3 Inclusion of First Order Corrections 



Here we examine what happens to the zeroth order classical solutions upon including the first 
order derivative corrections in the action. We already checked in [10] that the zeroth order 
Born-Infeld solution (C = 0) does not survive upon the inclusion of such corrections. More 
specifically, we numerically found a classical solution to the corrected field equations, but it 
was associated with an infinite value for the Lagrangian. Because as with zeroth order, the 
Lagrangian is 5*0(1, 1) invariant, the result that the of an infinite value for the Lagrangian 
follows for the entire orbit of solutions connected to the C = solution; i.e. case i). On the 
other hand, the case ii) BPS solution {\Q\ = \C\) is stable upon inclusion of the first order 
corrections, and just like at zeroth order, the Lagrangian vanishes. In fact the BPS solution is 
known to survive to all orders in the derivative expansion. [11] We shall verify that this result is 
consistent with the explicit expression for the first order terms obtained in [8], [9]. The stability 
analysis for the wormhole solutions case in) leads to the same results we obtain for case i). 
Namely, corrections to the zeroth order solution lead to an infinitely large correction of the 
Lagrangian. 

The first order corrections were initially computed in [8], [9] for the space-filling D9— brane. 

A dimensional reduction could then be performed to get the corrections to the DBI action 
(2.2) for an arbitrary Dp— brane. We first briefly recall the results of the dimensional reduction 
procedure at zeroth order. [4] One starts with the Born-Infeld (BI) action for the space- 
filling D9— brane. It is written in terms of a 10 x 10 matrix h with elements 

^AB = r/AB + 2W Fas , A^, - = 0, 1, ...,9 (3.1) 

where Fab = QaAb — QbAa is the ten dimensional field strength and we again assume a flat 
background metric tjab- h in (3.1) can be obtained from h in (2.1) by assuming the static 
gauge, which here means X"^ = for all A. S^j is given by 

4°] = - / d''^ 41 , 41 = 1 - \/-det[hAB] , (3.2) 

and from (2.3), Tg = l/{A7r'^a')^ . In dimensional reduction to the Dp— brane, 9 — p of the nine 
spatial directions are 'T-dualized'. Choose the T-dual directions tohe A = a = p+l,p+2, 9. 
One of the consequence of this procedure, is that the gauge potentials Aa in the T-dual 
directions get replaced with the transverse modes Xa of the Dp— brane according to 

2Tra'Aa ^ . (3.3) 
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The fundamental degrees of freedom are then X^^ and the remaining p + 1 gauge potentials 
Afj,., = 0, 1, which are functions of the p + I coordinates of the brane Then the 
nonvanishing matrix elements of h are h^^ and they are identical to h^j^i, of the DBI action 
written in the static gauge and given in (2.5). Finally after performing integrations in the 
T-dual directions (3.2) gets replaced by 



(3.4) 



where the Dp— brane tension is again given by (2.3), and one recovers (2.2) in the static gauge. 
So instead of working with h^^ as we did in the previous section we could have started with the 
10 X 10 matrix h. Then for the static spherically symmetric solutions of the previous section 
where just the p + 1 transverse mode is excited, 



gn 



i = 1,2,. ..p 



where f is the unit vector in the radial direction and spherical symmetry means / and g are 
only functions of the radial variable r. The 10 x 10 matrix h takes the form 

/-I -fr \ 
ff llpxp gr 
-gf 1 

V %-p)x(8-p)y 



h = 



(3.5) 



The first order corrections S^j to the action S^j of the space-filling D9— brane obtained 
in [8], [9] involve first and second derivatives of the field strength Fab- They are contained in 
the rank-4 tensor 



Sabcd = 2'Ka'dAdBFcD + {2TTa'fh^^{dAFcEdBFDG - dAFoEdBFcG) , (3.6) 
which is antisymmetric in the last two indices. Here h^^fiBC = Sq. The total action is 



o(0) , <^ 



(1) 

BI 



5/, 

gs J 



,10 



X < 1 



det[hAB] 



1 + -JA 



A = h^^hP'-'h^^h^'' {SbcegSdaij — 2SgibcSjeda) , 



(3.7) 



where n 



(27ra')S 
48 



. We again specialize to the case where a single transverse mode [the {p+l 



ith 



mode] is excited on a p < 8 brane, and consider static spherically symmetric fields. So the 
ansatz for h is again (3.5). Its determinant and inverse are given by 



det h 



-i + f-g^ 



(3.8) 



and 



h- 



det/t 



/I + 5' rf ~fg 

-ff -]l + (/2_^2)p 

-fg -rg f - 1 

V 



deth ]l(8_p)x(8-p)/ 



(3.9) 



13 



respectively. P is the projection matrix Pij = 5ij — fifj, satisfying Pijfj = and PijPjk = Pik- 
Some work shows that the nonvanishing components of Sabcd are 

Sijko = -Sijok = det h H/ hrjfk - ^ [llf + ^ {Pikfj + Pjkh) + Pijh 

Sijk p+i = -Sij p+i k = deth Hg' nrji-k - ^ (^g + ^ j {Pikfj + Pjkh) + j Pijh 

(Indet/i)' „ _ „..x 
^ijke = — ^ — [Pik^jve + PjeriVk - PurjVk - Pjknre) 

+ ^-^^^^^(PikPji - PiiPjk) , (3.10) 



where 



det h det h 

the prime here denoting derivatives in r. In addition we define 

Hk = fHg — gHf 

Substituting into the formula in (3.7) for A gives 



Hf = -^, Hg = -^ (3.11) 



+ ^ { (2 + (det hf^ {H/ - Hg' + Hk' 

+ - (In det hy' + - (In det hV - - (det hV 
2 r r 



+ ^{l + deth)(^p + l + {p-2)detf?j^ (3.12) 
So for the above ansatz the correction to the zeroth order Lagrangian density -C^^/ in (2.2) is 



'(1) 

(-det/i)3/2 



^'d'bi = . (3-13) 



In obtaining the equations of motion one must again write /(r) and g{r) in terms of 
potentials and cxtrcmize with respect to the latter. As the general system is quite involved, 
below we shall restrict to functions /(r) and g{r) which are related by a constant factor 

M . ^ , ,3.14) 

as what occurred for the zeroth order solutions (2.8). This set of configurations respects the 
50(1,1) symmetry (2.9) and (2.10). Using 

Ml^ap, «,„.o (3.15) 
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the Lagrangian density simplifies, and it is 5*0(1, 1) invariant. Once again tliere are three 
distinct orbits: i) \Q\ > \C\, ii) \Q\ = \C\ and Hi) \Q\ < |C|, and one expects that these orbits 
are classified by the corresponding value for the spatial integral of the Lagrangian density, 
which is now C,^dbi'^ ^^dbi- compute the latter we only have to examine one point on each 
of the orbits, which we do below. 

i) \Q\ > \C\. A convenient point on this orbit is the purely electrostatic case, where 
the transverse mode is suppressed: /(r) = 27ra'.4o(r) and g{r) = 0. Substituting this into 
^^DBi + '^DB/' ^^"^ varying with respect to Ao{r) gives the corrected Born-Infeld equation 



V— det h 



f 



det/i)3/2 



(-det/i)3/2 



2{p-l){3 + mrP-^Hfy 
(-det/i)3/2 



(3.16) 



+ 



(p - l)rP-^f 



f{p - 2 - r^i^l) + p{-2p - 3 + Ar'^Hj) + 3(-2p + 6 + ^r'^H}) 



(- det /i)3/2 

To get back the zeroth order equations set the left hand side equal to zero. So the right hand 
side represents the derivative corrections. The zeroth order Born-Infeld solution satisfies (3.16) 
as r — oo, so the corrections are negligible in this region. In [10], starting with the zeroth 
order Born-Infeld solution at r — >^ oo, we used (3.16) to numerically integrate to finite r. We 
found the resulting corrections to /(r) at finite r to be small, and just like at zeroth order, / 
tends to 1 as r ^ 0. Call fo{r) the zeroth solution, and /i(r) the correction it receives at first 
order. Below we plot fo{r) and /o(?") + fi{r) when Q = 1, C = and p = 3: 




fig.7 /o and /o + /i vs r 

f or p = 3, Q = 1, C = 

In the above we set k = 1 which is equivalent to choosing the scale for r. If the zeroth 
order Born-Infeld solution gives a reasonable approximation to a classical solution in the full 
effective theory, and one can apply the derivative expansion to get the latter, then higher order 
corrections should be small. In particular, we expect only a small change in the value of the 
Lagrangian at the next order. If one assumes this to be the case a Taylor expansion about the 
zeroth order solution gives 



15 



h (3.17) 

/=/() 



5/ 

The last term vanishes by the field equations, and so the first order correction to the Lagrangian 
is J [-^Ds/K/o)- However, it is easy to check that ['C^^j-](/o) diverges near the origin as 
I/t-Sp+i gQ ^jjg gj,g^ order correction is not small; Rather, J dP^ [i2^^j](/o) is infinite! This 
agrees with the result found in [10], and indicates that the Born-Infeld solution, and indeed all 
case i) solutions, are unstable under inclusion of first order derivative corrections. 

a) \Q\ = \C\. This is the BPS case g{r) = f{r). Here -C^^j = 0, and so just like at 
zeroth order the Lagrangian vanishes. To find equations of motion we must first vary / and g 
(or more precisely Xp+i and ^0) separately and then impose the BPS condition. We do this 
in Appendix B. (Actually, there we don't impose the restriction of spherical symmetry.) The 
result is simply 

|i + 2k (V2)4 ^0 = , (3.18) 



with the same equation for Xp+i. For the case of spherically symmetric solutions we can use 
= :p^drrP~^dr. Eq. (3.18) says that the zeroth order solution is also valid at first order. 
This agrees with [11], where it was shown that the BPS solution is valid to all orders. The 
result (3.18) thus provides a check of the computations in [8], [9]. 

in) \Q\ < \C\. A convenient point is the purely transverse case. Here the electric field 
vanishes: /(r) = and g{r) = X^^^(r). Substituting this in (3.13) and varying Xp+i(r) gives 



■\/— det h 



rP~^g - C\J- det ~h)j 



(-det ^)3/2 



(-det/i)3/2 (-det^)3/2 ^ ■ ' 



det ft.)3/2 

Again to get back the zeroth order equations set the left hand side equal to zero, and so 
the right hand side represents the derivative corrections. (3.19) is also obtained by making 
the transformation /(r) igir) and Q — > iC in (3.16), so solving for real g(r) in (3.19) is 
equivalent to solving for imaginary /(r) in (3.16). As with case i), starting with the zeroth 
order solution at r ^ 00, we can use (3.19) to numerically integrate to finite r. We find that 
just like at zeroth order, g becomes singular at some finite r, which appears to be slightly 
greater than tq. We can then conclude that the corrections cause the wormhole to become 
wider. Call go{r) the zeroth order solution, and gi{r) the correction it receives at first order. 
Below we plot go{r) and go{r) + gi{r) when C = 1, Q = and p = 3: 
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1.5 



2.5 



fig.8 go and go + gi vs r 

for p = 3, C= 1, Q = 

Again we set «; = 1. Figure 8 shows that the correction gi to the solution is small away from 
the wormhole throat. On the other hand, the corresponding correction to the Lagrangian 



density appears not to be small, as is indicated below where we numerically compare [C 



^DBiKao+gi) and [/:^^j](5o): 



(0) 



(0) 

DBI 



+ 




fig.9 



i^DBiKdo) and [i:^^^ + C)^'sj]{go + gi) vs r 
f or p = 3, C = 1, Q = 

^(1) 



,(0) 



The lower curve in figure 9 is [-C^^^ + jC^j^'^jjigo + 9i)- We then expect large corrections to its 



integral. In fact, we find that the numerical integration of [-C^^^ + i2^^j](5fo + fails to give 
a convergent result. So just as in case i), the integral of the Lagrangian density appears to 
be ill-defined, indicating that the case in) solutions are also unstable under inclusion of first 
order derivative corrections. 



4 Conclusion 

The preliminary indications here are that the classical wormhole solution may not be a reason- 
able approximation to a solution in the full D— brane theory. If so it therefore cannot prevent 
the decay of the brane-anti-brane system. More generally, it appears that the only solution 
that survives higher order derivative corrections may be the BPS solution. On the other hand, 
a more extensive analysis may be possible. For example, it would be interesting to drop the 
restriction to static solutions. Perhaps time dependent configurations can survive first order 
derivative corrections, or perhaps the zeroth order static solutions evolve to time dependent 
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ones after including the higher order. To check this would require combining two separate sta- 
bility analyses, which we referred to as a) and b) in the introduction. A final but unpleasant 
(from a computational point of view) possibility is that solutions are recovered only after going 
beyond the first order. Moreover, all orders may be required, meaning the solutions may be 
non perturbative. 
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Appendix A Alternative Gauge 

Here we reconstruct the zeroth order wormhole solutions starting with the action written in 
an alternative gauge. This gauge is obtained by replacing the r— coordinate of the static gauge 
by 2; = Xp^i gauge. It has the advantage that it removes the coordinate singularity appearing 
at the midway point of the wormhole, and shows that there is a smooth solution connecting 
the brane to anti-brane. Coordinate singularities re-appear, however, at the location of the 
brane and anti-brane. We show that the electric field is well behaved in this gauge, and in fact 
is a constant. 

Consider the domain S^^^ x M?, with local coordinate patches used to parametrize the 
Dp— brane. Denote by z one of the coordinates of M?, while the other corresponds to time t. 
We look for spherically symmetric static solutions with 

Xo = t, Xf + Xi + --- + X^ = R{zf, Xp+i = z (A.l) 

So for example for the D3— brane we can write 

Xi = R{z) sin 6 cos (j) 
X2 = R{z) sin 6 sin (j) 

X3 = R{z) cos e , (A.2) 

using standard spherical coordinates and <j), 

htt = -l, hee = R{zf, h^^ = R{zfWe, K^ = R'{zf + l, (A.3) 

where here the prime denotes a derivative in z. Now introduce a 2;— dependent electrostatic 
field in the z— direction, leading to the off-diagonal components 

htz = -h,t = -E{z) (A.4) 

In terms of the electrostatic potential ^0 which we now write as a function of z, E[z) = 
2iTa' dzAQ{z). After performing the angular integrations, the DBI Lagrangian L^j^^j (ignor- 
ing the vacuum term) will be proportional to R{zY \/R'{zY — E{zY + 1. Generalizing to 
arbitrary p, 

oc Riz)P-^ ^/R'{zf - E{zf + 1 (A.5) 
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Variations in the electrostatic potential and R{z) give 

R{z)"-' E{z) 



R"{z)R{z) = {p-l) {Rl{zf - E{zf + 1) 
respectively. From the first equation 

E{z) 



Q 



' R'{zY + 1 



(A.6) 
(A.7) 

(A.8) 



where Q is an integration constant, and substituting into the second equation 



After integrating once 



1 



(A.9) 



(A.IO) 



where C is an integration constant. Since R'{z) corresponds to l/g(r), the result agrees with 
(2.8). For the wormhole solution, R{z) is nonsingular everywhere except at the location of the 
branc and anti-brane. At the midway-point on the wormhole, i? is a well-defined function of 
z, and is a minimum since 



(A.11) 



and \Q\ < \C\ for wormhole solutions. So now coordinate singularities appear at the brane and 
anti-brane, rather than at the midway-point on the wormhole. By substituting (A.IO) into 
(A.8) one gets that the electric field E{z) in the z— direction is a constant 



\E(z) 



(A.12) 



which agrees with (2.23). It goes to one in the BPS limit, and l/\/p — 1 for the minimum 
energy wormhole. We conclude that in this coordinate frame there are no singularities in the 
electric field (for Cj^O). 



Appendix B 



First order BPS equations 



Here we derive the first order BPS equation (3.18). Unlike in section 3, we make no 
restriction to spherical symmetry. Our starting point is then not (3.5), but 



f-l -f 
f ^pxp 9 
-9 1 

V 



ll(8_p)x(8-p)/ 



(B.l) 
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(B.2) 



where / = 2'iTa'VAo 9 = ^^p+i are vector fields on the Dp— brane. The general BPS 
condition is f = g. Upon imposing this condition on one gets 

/-l-P -f f\ \ 

h \bps=\ ^ ^ 

\ %-p)x(8-p)y 

while the only nonvanishing components of Sabcd are 

Sijko\BPS = SijkpJr\\BPS = —SijOk\BPS = —Sijp+l k\BPS = didjfk (B-3) 

Since the BPS action vanishes, we must impose f = g after performing the variations in the 
action to find the BPS field equations, i.e. we must be allowed to perform separate variations 
of ^0 and Xp+i . By varying and then setting f = g , 





BPS = 


-2 + 




BPS = 


-2 P /• df 




BPS = 


{l + 2p)f.5f = <5^^+^%P5 




BPS = 


(1 + f ) 8fi + / ■ 5ffi = -Sh'^'lBPs 




BPS = 


-P Sfi - / ■ Sffi = -dhP+' %ps 




BPS = 





(B.4) 

while the nonvanishing components of 6 Sabcd are 

SSijkolBPS = —SSijoklBPS = didjSfk + SSijkp+i\BPS 

SSijkp+i\BPS = -SSijp+i k\BPS = difi djfk Sfi - hP^diSfk djfi + {i ^ j) 

SSijkelBPS = djfi di6fk + difk djdfe - {i ^ j) (B.5) 
In evaluating 6A\bps we can use (h^^ SijAB)\BPS = 0- Then 

5A\bps = -^hP^'Sj^DA {5h^''Si,BC + h^^'^SijBC + h^^Sh'^SejBc} \bps 

= -SdidjfkdidjSfk (B.6) 

Now substitute f = 27ra'Vylo to obtain (3.18) from the variation of Ao in jC^^^^j ''r ^d^bi- 
gets the same results from variations of Xp^i . 
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